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INTRODUCTION 


1. 0 

The  purposes  of  this  work  are  to  provide  an  analytical  basis  for 
determining  the  effects  of  phase  distortion  in  linear  bandpass  systems  and  to 
develop  a model  of  a class  of  linear  bandpass  system  functions  which  is  use- 
ful in  numerically  estimating  system  phase  0(w),  and  its  derivative,  group 
delaY  [tgrM]- 

The  motivation  for  this  work  stems  from  consideration  of  the 
Aerosat  wideband  channel  group  delay  specifications  in  the  period  November  - 
December,  1975.  Two  applications  notes  pertinent  to  that  system  are  given 
in  Sections  4.  1 and  4.  2. 

Although  the  results  of  this  work  are  generally  applicable  to  lowpass 
and  highpass  system  functions,  primary  interpretation  is  carried  out  in  terms 
of  the  bandpass  function  shown  in  Figure  1.  0-1. 

In  Section  2.0,  paried-echo  theory  is  developed.  This  provides  a 
superposition-of-modulating-signal  technique  for  determining  signal  distor- 
tion once  the  phase  nonlinearities  are  known.  In  Section  3.  0,  the  phase  non- 
linearities  for  a minimum -phase  transfer  function  of  the  type  shown  in  Fig- 
ure 1.0-1  are  determined. 

Section  5.  0 is  a discussion  and  comparison  with  previous  work  and 
Section  6.  0 restates  the  principal  conclusions. 

2.  0 PAIRED-ECHO  THEORY 

The  purpose  of  this  section  is  to  provide  an  analytical  basis  for 
answering  the  questions: 

• How  can  the  degradation  of  a signal  due  to  nonlinear 
phase  be  computed? 

• What  is  the  basis  for  AM  to  PM  and  PM  to  AM 
conversion  ? 
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2.  1 The  Theori 


The  system  function  indicated  in  Figure  1.  0-1  is  modeled  as  shown 
in  Figure  2.  1-1;  the  usual  relationships* 


6(w) 

:"J  |H(w)|, 


= H(w)  X(w), 
= h(t)  x(t) 


= convolution  operator). 


hold. 


The  phase,  G(uj),  is  expanded  in  a Taylor  series  about  the  center 


frequency,  wQ: 


e(u>)  - 0Q  + ( ui-  ) 1 0 + 0 (w  — u)Q ) , 


t ' t " t 

e(w  - uio)  = -yr  (w  - wQ)2  + “fy  (w  - wQ)3  + • • • + ° n!  (w  - w0)n» 


where 


0Q  = 0(wQ)  = carrier  delay  at  center  frequency, 

tQ  = ^ 4 r(w0)  = group  delay  at  center  frequency,  and 

o 

( ) d^n) 

t ' = : — r t (co)  = derivatives  of  the  group  delay  at  center  frequency. 

dur*  ^r 


* _-0(w)  :6(w) 

The  convention  H(w)  = |H(w)je  J as  opposed  to  H(w)  = |H(w)|eJ  is 
chosen  for  convenience  in  dealing  with  the  (negative  phase)  bandpass  system 
function. 


In  Appendix  1,  it  is  shown  that,  for  a 0(w  - wQ)  representable  as  a Fourier 
series  in  the  interval  (w^,  w3)  (only)*,  the  output  of  the  system  shown  in  Fig- 
ure 2.  1-1  when  subjected  to  an  input 

;wt 

x(t)  = a(t)  eJ  , a(t)  real  or  complex  (2.  1-4) 

is 

y(t)  = ye(t)  ® O'1  (|H(oo)|U(w)],  (2.1-5) 

where 

= inverse  Fourier  transform  operator  and  U(w)  is  as  shown 
in  Figure  2.  1-2. 

The  part  |H(w)|U(w)  of  (2.  1-5)  is  the  modification  of  the  system  out- 
put due  to  the  magnitude  response  of  the  system,  H(w)  , and  the  interval 
(w2,  w ) over  which  the  Fourier  expansion  for  0(oj  - gjq)  is  valid.  In  most  of 
what  follows,  attention  is  restricted  to  y0(t).  Any  degradation  which  occurs 
in  y ( t ) is,  in  baseband  terms,  simply  lowpass  filtered  by  the  lowpass  ver- 
sion  of  iH(w)|U(w). 

In  Appendix  1,  it  is  shown  that 


V‘>  = 


» :c  \ .[wt-e0+(w0-w)(t0+nT)] 

= 2^  CneJ  - ‘o  " nr)’ 


(2.  i-7; 


n=-a> 


where 


2rr 


’3  " “2  3 *2 


-,  and 


‘^3  _.[e(w-wo)+2Trn(u)-wo)(u)3-w2)J 


(2.  1-8) 


n u3  - 


du>. 


The  Taylor  the  i Fourier  series  expansion  of  0(w  - wQ)  proves  convenient 
when  the  properties  of  the  Fourier  coefficients  are  considered  (Section  2.  3). 


To  the  extent  that  the  interval  is  finite,  the  result  (2.1-5)  is  approx- 

imate; in  practice,  it  can  occur  that  tractable  expansions  for  H (w)  are  only  known 
in  some  finite  bandpass  interval  (w^.oj^).  Then,  (2.  1-5)  is  only  that  part  of 

the  output  which  is  due  to  the  bandpass  interval.  To  this  must  be  added 
that  part  of  the  output  due  to  skirt  response.  As  will  be  exemplified  in  Sec- 
tion 3.  0,  it  is  oftentimes  convenient  to  separate  in-band  and  out-of-band 
responses. 


Heuristic  Interpretation 


For  the  purposes  of  initial  interpretation,  it  is  assumed  that  w = u»Q 
and  that  a(t)  and  the  cn  are  real  in  (2.  1-9).  The  case  of  a complex  a(t)  is 
treated  later  in  the  present  section,  and  the  case  of  complex  cn  is  treated  in 
Section  2.  4. 


In  order  to  illustrate  the  effects  of  phase  distortion,  consider  the 
following  example:  Let  a(t)  = ±1  (biphase  modulation)  with  period  3t  and  let 
the  phase  distortion,  0(w  - w ) be  such  that  c =0.  75,  c.  = 0.  25,  c , = 0.  25 

r o'  o 1 - 1 

cn  = 0,  n 4 0,  1,  -1..  The  modulation  present  at  the  output  of  the  phase  box 

is  shown  in  the  lowest  diagram  of  Figure  2.  2-1.  Clearly,  this  is  not  a 

replica  of  that  at  the  input.  Depending  on  the  specific  nature  of  the  detector 

this  could  cause  bit  errors  at  the  receiver. 


If  there  was  no  phase  distortion,  i.  e.  , if  6(w  - wQ)  = 0 in 
(2.  i-2),  c = 1,  c = 0,  n 4 0.  Then,  the  output  modulation  is  a delayed 
(by  t ) version  of  the  input  modulation,  a(t)  (second  diagram  in  Figure  2.  2-1 
with  c = 1 4 0.75).  In  the  general  case,  however,  there  are  an  infinite  num- 
ber of  superposed  signals  due  to  phase  nonlinearity. 

The  above  interpretation  is  easily  extended  to  PM  signals.  Let 

•4>(t) 

a(t)  = eJ  . (2.2-1) 


Then,  yQ(t)  is,  from  (2.  1-9), 


Vt> 


t-eo) 


.4>(t-t  -nr) 

V cn  e->  ° 


(2.  2-2) 


The  output  modulation  (following  the  2 sign)  is  still  delayed  by  tQ  + nr;  in 
some  works,  group  delay  is  called  "envelope  delay.  " In  the  present  example, 
there  is  no  varying  envelope  but,  yet,  the  modulation  is  delayed  by  tQ  + ht. 
For  obvious  reasons,  the  term  "group  delay"  is  to  be  preferred. 

The  graphical  interpretation  exemplified  by  Figure  2.2-1  is  still 
applicable.  In  the  PM  case,  one  must  consider  both  the  real  and  imaginary 
parts  of  the  output  modulation. 

2.  3 Properties  of  the  cn 

By  the  change  of  variable 


4 = 


(2.3-1) 


in  (2.  1-8), 


9 


_1_  r .j[e(K|).n5] 


n 


JL 

2tt 


2n 


d£, 


K = 


“3  ' “2 
2it 


(2.  3-2) 


cos  [e(K£)  - n£]  d|  - j / sin  [e(K£)  + n|]  d£ 


L^-it 


'-IT 


where,  from  (2.  1-2), 


e(K£)  =^-(K£)2+-^-(K£)3+  ...  +-^  (K£)n.  (2.3-3) 


t ln~l) 


In  Appendix  2,  it  is  shown  that, for  any  minimum-phase  system  which  is  sym- 
metric in  magnitude  about  with  a bandwidth  much  less  than  twice  the  center 

frequency,  t (u>)  is  even  about  i.  e.,  if 
8^  * 


|H(Wq  + Aw)  | = | H(w0  - Aw)i 

and  2wq  » - w-,, 


(2,  3-4) 


then 


tgr  (wQ  + Aw)  = tRr  (u>Q  - Aw). 


(2.3-5) 


It  follows  from  this  result  that  the  odd  derivatives  of  tQ  = tgr(w0)  are  zero: 


t^  = 0,  n odd. 


(2.  3-6) 
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(2.  3-7) 


Then,  from  (2.  3-3)  with  tQ^n^  = 0,  n odd, 

0(-K£)  = -e(K£). 


That  is,  9(K£)  is  odd. 
assumes  the  form: 


-j 


The  second  term  of  the  second  expression  in  (2.  3-2) 


sin  [odd  function  of  d£  = 0. 


(2.  3-8) 


Then,  the  cn  are  real.  By  expanding  cos  [e(K4)  - n£]  as  the  sum  of  sin  sin 
and  cos  cos  terms  and  making  use  of  the  property  (2.  3-7),  it  is  easy  to  show 
that,  for  the  restrictions  (2.  3-4), 


c 

-n 


"Cn  ~ 2tt  / sin  9(K£)  sin  n£  d4>  real.  n^O. 


(2.  3-9) 


Now,  attention  is  directed  to  another  property  of  the  cn:  convergence. 
Causality  implies  that  the  system  cannot  respond  before  an  input  is  applied. 
In  terms  of  the  example  in  Figure  2.  2-1,  the  n in  the  c is  bounded  to  a 
maximum  of  2;  otherwise,  the  output  would  occur  before  the  input.  In 
general, 


c = o for 
-n 


n > lowest  integral  value  of 


(2.  3-10) 


At  least  in  the  case  (2.  3-4),  this  also  bounds  the  maximum  n in  c via  (2.  3-9). 

n 


11 


- HUM 


..  ...  . 





« ....... 


; 


Finally,  attention  is  directed  to  the  uniqueness  of  the  cn:  for  each 
choice  of  the  interval  (w^,  w^).  there  exists  a different  set  of  cn  spaced  at 
t = 2tt/(u>2  - oj^)  seconds.  The  choice  of  t is  somewhat  judgemental.  Often- 
times, there  is  a tractable  expansion  forH(w  - wQ)  applicable  to  the  bandpass 
interval  only  (e.  g. , the  3dB  or  ldB  bandwidth)  which  will  dictate  the  choice. 

It  should  be  noted  that,  as  the  interval  t becomes  larger,  a larger  set  of  the 
cn  have  to  be  computed  [See  (2.  3-10)  and  (2.3-9)]. 

2.  4 A Comment  on  AM  to  PM  and  PM  to  AM  Conversion 

In  the  general  (asymmetric  |H(uj)|)  case,  the  cn  are  complex.  From 
a term-by-term  expansion  of  the  output  (2.  1-9)  of  the  phase  box  involving 
Re  c + ilm  c , it  is  easy  to  show  that  PM  is  induced  in  an  AM  signed 
[a(t)  real]  and  that  AM  is  induced  in  a PM  signal  [a(t)  = eJ®'  '].  Both  effects 
are  undesirable.  Therefore,  arithmetically  symmetric  systems  [cn  realj 
are  desirable.  • 

By  using  (2.  1-7)  and  expanding  in  the  manner  indicated  above,  it  is 
easy  to  show  that,  even  if  |H(uj)|  is  symmetric,  the  aforementioned  AM  to 
PM  and  PM  to  AM  conversion  effects  occur  when  the  signal  center  frequency, 
u>,  is  not  equal  to  the  system  center  frequency,  u)Q. 

3.  0 GROUP  DELAY  OF  AN  IDEALIZED  MINIMUM-PHASE 

BANDPASS  SYSTEM 

In  the  writing  of  system  specifications,  the  magnitude  character- 
istics are  usually  considered  first.  These  control  out-of-band  rejection, 
smoothness  of  the  in-band  response,  intermodulation  distortion,  etc.  If  only 
voice  communications  with  standard  non-encoded  AM  and/or  PM  modulation 
is  employed,  amplitude  specifications  are  all  that  is  necessary,  as  the 
human  ear  is  substantially  oblivious  to  phase  distortion.  However,  if  digital 
and/or  navigation  signals  are  employed,  group  delay  is  usually  considered. 
Oftentimes,  a tradeoff  must  be  made  between  phase  and  amplitude  specifica- 
tions and  the  cost  of  equalizing  (i.  e.,  making  more  linear)  the  phase 
characteristics. 


! 
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The  thrust  of  the  present  work  is  that  of  taking  as  given  the 
amplitude  characteristic  and  inferring  the  phase  characteristic.  This 
addresses  the  practical  questions  of,  e.g., 

• How  does  (amplitude)  skirt  steepness  affect  group 
delay? 

• How  does  in-band  (amplitude)  ripple  affect  group 
delay? 

In  the  following,  minimum-phase  (no  zeros  in  the  right  half  plane) 
filters  are  considered.  In  practical  terms,  this  means  that  unequalized 
filters  are  considered.  This  includes  Butterworth,  Tchebyshev,  Elliptic 
function  and  other  such  filters  that  are  normally  used  to  meet  amplitude 
specifications.  Even  if  it  is  intended  to  equalize,  this  minimum-phase  prob- 
lem is  generally  considered  first  in  order  to  assess  the  difficulty  of  design- 
ing the  equalizer;  most  "equalized”  systems  are  a cascade  of  a minimum- 
phase  network  and  an  "all-pass"  equalizer. 

Again,  Figure  1.0-1  indicates  the  class  of  systems  under  considera- 
tion. In  Section  3.  1,  the  case  of  a flat  in-band  magnitude  response  is  treated. 
Section  3.  2 extends  this  to  an  arbitrary  magnitude  response, 

3.  1 Group  Delay  of  an  Idealized  Minimum-Phase  System 

Without  In-Band  Amplitude  Distortion 

The  principal  aim  here  is  to  restate  the  results  of  Holmes  and 
Keenan^  ^ pertinent  to  the  filter  shown  in  Figure  1.0-1  with  a flat  in-band 
response.  These  will  be  used  later  in  the  applications  notes  in  Section  4.  0. 

Holmes  and  Keenan  show  that,  for  the  above  bandpass  system, 


t (w)  = - 
gr'  ' it 


In 


"|H(cj2)|" 

i 

|h(Wi)| 

w 2 “ W1 

_ 

_ 

U) 

(3. 1-1) 
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K , 


l 

I 4 


1 


1 


It  is  also  shown  that,  if  2wQ  » (u^  - w^), 


t (w)  = 
gr'w; 


+ — L-ln|i+  4 


J4  w3 


^2  ^ w ^ ^3- 


In  the  symmetric  case  |H(w  + Au>)|  = [H(w  - Aw)[,  (3.  1-2)  reduces  to 


g 


r<“>  “ ^ln  [|h|u^|]  ln  [( 


1 + 


ft 


1 + 


ft 


L_) 

- Co/  , 


OJ  ^ U)  ^ • 


- u2  " wl  = w4  ” w3’ 

Now,  one  measure  of  the  amount  of  group  delay  distortion  is 


From  (3.  1-3),  this  is  easily  shown  to  be 


^ t 


|H(u,2)|l 

(l  + 1 + k)  (l  + 1 - k) 

jH(Ul)|Jl 

d + a)2 

where 


(3.  1-2) 


(3.  1-3) 


At^r  = maximum  group  delay  - minimum  group  delay  (3.  1-4) 


(3.  1-5) 


a = 2n/(u>2  - ^2)  > 0 and 

K = Fraction  of  filter  bandwidth  occupied  by  the  signal  < 1. 


(3.  1-6) 


The  case  K = 1 corresponds  to  the  signal  bandwidth  equal  to  the 
filter  bandwidth  and,  due  to  the  discontinuous  slope  of  the  magnitude  charac- 
teristic occurring  when  the  in-band  response  is  flat  in  lo^)  in  Figure  1.0-1, 


l 

4 


11+ 


inlinite  At  . In  a practical  approximation  of  the  magnitude  characteristic 
6 

(Butterworth,  Chebyshev,  etc.  ),  Atgr  is  not  infinite.  The  model  tends  to 
provide  estimates  of  At^r  which  are  large  in  comparison  with  those  incurred 
in  a practical  realization;  the  closeness  of  the  upper  bounding  decreases  with 
increasing  K. 

3.  2 The  Effects  of  In-Band  Amplitude  Distortion 

In  this  section,  the  work  of  Holmes  and  Keenan  is  extended  to 
include  the  effects  of  an  arbitrary  continuous  in-band  amplitude  function, 

F(oj  - ^2).  In  Appendix  3,  it  is  shown  that  this  function  can  be  represented 


E(w  - Wz)  = In  ^ 


CO 

■E 


b sin  nTT 


\w3  " w2 ) 


(3.2-1) 


where 


n u>2  - u»2 


F(ui  - u ) sin  nn 


/ w"“2  \ 

\w3  " w2/ 


(3.  2-la) 


Figure  3.2-1  illustrates  some  special  cases;  all  b^  = o,  n even,  corresponds 
to  the  symmetric  |H(uj)|  case. 

In  Appendix  2,  it  is  shown  that  that  part  of  the  group  delay  due  to 
(arbitrary)  amplitude  distortion  of  the  form  (3.  2-1)  is 


F(u>  - u >2)  beings  expanded  about  lo 2 in  terms  of  a sin  function  only  does  not 
restrict  the  eveness  or  oddness  about  u>0,  but  only  the  continuity  of  F(u>2): 
F(0)  = F(u)3).  (3.  2-1)  is,  then,  an  arbitrary  but  continuous  function. 
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=£  2 
, n ¥ 3 

: o,  n 4 1 


1 Cases  of  Amplitude  Distortion:  Ripple 
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nuAu) 


niT  Aoj 


a - bAw)+  Si  (a  + 


where  the  only  approximation  which  has  been  made  is  2<jjq  » and 


In  the  special  case  of  arithmetic  symmetry,  the  terms  following  the  second 
summation  in  (3.  2-2)  are  zero.  In  the  special  case  of  the  ripple  indicated  in 
Figure  3.2-1,  all  terms  but  one  in  (3.2-2)  are  zero. 

The  group  delay  incurred  in  the  interval  (u>2>  w^)  is  that  in  (3.  2-2) 
plus  that  computed  in  (3.  1-1). 


4.  0 APPLICATIONS  NOTES 

Applications  of  the  foregoing  theory  were  incurred  in  the  course  of 
specifying  the  (wideband)  Aerosat  satellite  transponder.  For  a narrowband 
application,  see  Holmes  and  Keenan.  ^ 

4.  1 Differential  Group  Delay  of  the  Aerosat  Wideband  Channel 

In  writing  a specification  for  the  Aerosat  wideband  (=  10  MHz)  chan- 
nel, it  was  desired  to  specify  the  spacecraft  differential  group  delay  in  such 
a manner  as  to  not  force  the  task  of  equalization  on  the  spacecraft  designer; 
equalization  (if  nece ssary)  is  best  carried  out  at  the  ground  receiver  sites. 
The  problem,  then,  is  to  determine  the  group  delay  incurred  if  a minimum- 
phase  (i.e.  , unequalized)  filter  is  used  to  meet  the  satellite  magnitude 
response  requirements. 

The  sole  constraints  on  the  magnitude  response  of  the  satellite  were 
that  the  noise  bandwidth  was,  maximally,  10  to  12  MHz  and  that  a signal  of 
bandwidth  between  5 and  8 MHz  had  to  be  handled  with  little  or  no  amplitude 
degradation. 

The  channel  characteristics  were,  the,  of  the  nature  shown  in  Figure 
1.  0-1  with  arithmetic  symmetry  and  a flat  in-band  response.  In  Appendix  4, 
it  is  shown  that  the  noise  bandwidth  of  this  channel  is  * 


N a <f3  - fz> + hrrfoooi 


n = f„ 


fl  = f4  - V 


(4.  1-1) 


With  N fixed  at  10  MHz  or  12  MHz,  a relationship  is  established 
between  U and  the  bandwidth,  (f3  - f^).  This  relationship  is  used  in  (3.  1-5) 
to  obtain  the  results  given  in  Figure  4.  1-1. 


* , 

The  approximation  arises  due  to  neglecting  the  contribution  to  the  noise 
bandwidth  outside  the  60dB  bandwidth  of  the  filter. 


NOISE  BANDWIDTH 

10  MHz 

12  MHz 


SIGNAL 

BANDWIDTH  (MHz) 


BANDWIDTH  OF  FILTER  (f MHz 


Differential  Group  Delay  Incurred  in  an  Idealized 
Unequalized  Channel  with  Fixed  Noise  Bandwidth 
as  a Function  of  the  Channel  Bandwidth 


For  the  5 MHz  signal  bandwidth  case,  it  was  recommended  that 

30  nsec  At  be  allowed  and,  for  the  8 MHz  signal  bandwidth  case,  that  a 
£>  ^ 

At  of  150  nsec  be  allowed.  This  permits  the  spacecraft  designer  to  use 
8 * 

unequalized  (minimum-phase)  filters  to  meet  the  more  desirable  noise  band- 
width specification  of  10  (as  opposed  to  12)  MHz. 

4.  2 Degradation  of  Autocorrelation  Signals  Due  to 
Phase  Distortion 

In  a class  of  systems  of  practical  interest  (e.  g. , Aerosat),  the  out- 
put of  a linear  channel  is  multiplied  by  a locally-generated  replica  of  the 
modulation  impressed  on  the  carrier  in  order  to  recover  data.  If  the  modu- 
lation on  the  carrier  is  a(t)  (where  a(t)  is,  in  general,  complex),  the  base- 
band signal  is,  for  a linear-phase  system  with  infinite- bandwidth  magnitude 

j 

response 


0(t  - t0  - M + n(t)]  [a(t  - a)]  = a(t  - tQ  - \)  a(t  - a)  + n(t)  a(t  - c r), 


where 


(4.  2-1) 


n(t)  = noise, 

tQ  = group  delay  = constant,  and 
X.  = unknown  transit  time  of  the  signal. 


(4.  2-2) 


and,  in  practice,  a is  varied  (usually  in  a delay-locked  loop^)to 


maximize 


the  function 


R (<r 
Y 


to  " “ 2T 


r 

Ljl 


a(t  - t - X.)  a(t  - tr)  dt, 


(4.  2-3) 


Modification  of  a(t)  by  the  finite  bandwidth  of  the  system  is  not  under  con- 
sideration here;  only  the  degrading  effects  of  phase  nonlinearity  are  being 
considered.  lH(u>)|  = 1 for  this  purpose. 


20 


where  2T  is  the  period  over  which  integration  (integrate  and  dump)  takes 
place.  The  maximum,  Ry(0),  occurs  when  <r  = tQ  + X.  Because  tQ  is  known 
a priori,  determination  of  <r  provides  a measure  of  X,  hence  the  transit  time 
of  the  signal  and  navigation  and  other  data.  If  the  integration  period,  2T, 
extends  over  many  intervals  of  t(t  = 1 /system  bandwidth),  R^fcr  - tQ  - X)  is  a 
good  approximation  of  the  formal  autocorrelation  function  of  a(t  - t - X). 


The  question  arises  as  to  the  effect  of  nonlinear  phase  on  the  auto- 
correlation function  of  the  output  of  the  system.  Consider  the  signal  part  of 
the  output  translated  to  baseband  [w  = 0 = 0 in  (2.  1-9)]: 


For  the  moment,  assume  that  the  system  is  minimum-phase  and 
symmetric  in  magnitude  so  that  c n = -cn,  real  [See  (2.  3-9)J.  Then,  for  a 
nonlinear  phase  system,  the  autocorrelation  function  of  the  output  of  the  sys 
tern  (i.  e.,  the  output  "signal")  is 


where  the  integrals  following  the  (one-sided)  £ stem  from  c 


The  autocorrelation  of  the  input,  a(t  - X),  is  defined  as 


R (<r  - X) 
x 


a(t  - X)  a(t  - X - cr)  dt. 


(4.  2-6) 


By  comparison  of  the  first  term  of  (4.  2-5)  with  (4.  2-6),  it  is  appar- 
ent that  the  first  term  of  (4.  2-5)  can  be  written  as 


c R (<r  - X - t ). 
o x o' 


(4.2-7) 


With  the  changes  of  variable  = t - tQ  - nT,  ^ = * " nT  *n 

first  and  second  integrals  enclosed  in  the  brackets  of  (4.2-5),  respectively, 
these  terms  can  be  written  as 


T-t  -nT 
o 


/-T  -t  -nT 
o 


T -t  +nT 
o 


- X)  a(4j  - X - 0-)  d^  - 


■ T-t  +nT 
o 


a(| 2 “ M a(^2  ~ ^ ^2’ 

(4.2-8) 


If 


T » t + nT, 
o 


(4.2-9) 


(4.2-8)  is  approximately  zero. 

From  (2.3-10),  it  is  known  that  c|n|  = 0 for  n > tQ/T  (where  the  Ini 
in  cn  follows  from  (2.  3-9)).  Thus,  n < tQ/r  (i.  e.  , nr  < tQ)  as  regards  com- 
putation of  the  terms  (4.  2-8).  The  condition  (4.  2-9)  can  be  written  as 

T » t + t = 2t  . (4.  2-10) 

o o o 

In  most  practical  systems  (including  Aerosat),  this  condition  is  met. 
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If  the  condition  (4.2-10)  prevails,  it  follows  from  the  above  that 

Ry(<r  - \ - tQ)  s cQ  Rx  (cr  - \ - tQ)  (4.  2-11) 


That  is  to  say  that,  if 


a)  The  system  is  minimum-phase  and  symmetric  in 
magnitude  about  , (|H(  wq  + Aw ) j = | H(  w - A w)  | ) and 

b)  The  input  signal  center  frequency  is  equal  to  the 
system  center  frequency  (w  = u)Q)  and 

c)  Twice  the  center  frequency  is  much  greater  than 


the  bandwidth  [2w  » ( 


w3  " “2 


)]  and 


d)  The  integration  time,  2T,  is  much  greater  than 
four  times  the  midband  group  delay, 

4t0  (2T  » 4t0>, 


(4.  2-12) 


the  sole  degradation  of  the  autocorrelation  signal  due  to  nonlinear  phase 
(only;  recall  that  the  effects  of  the  magnitude  response  have  not  been  con- 
sidered) is  expressible  as  (4.2-1).  The  constant  cq  can  be  computed  from 
(2.3-2): 


2tt 


cos 


J.II  (llll 

(K 4 )3  +~  (K4)5  + ...  + 


t (h-1) 


n: 


(K|) 


n 


di  < i. 


K = 


w3  - w2 


- f - f 
2 it  *3  12* 


The  above  is  not  to  say  that  one  need  not  consider  the  effects  of  non- 
linear phase  on  autocorrelation  signals,  but  only  that  these  should  be  con- 
sidered from  a different  viewpoint  than  that  which  is  conventionally  employed. 
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First,  arithmetic  symmetry  of  the  magnitude  characteristic  should  be 

/ \ 

insured.  One  (of  many)  was  of  doing  this  is  to  bound  t ' , n odd. 

Then,  from  (3.6-10),  tQ^n\  n even,  should  be  bounded  to  minimize 
the  signal  loss  (i.e.,  maximize  cQ).  "Differential  group  delay" 
becomes  an  obscure  measure  of  performance. 

If  one  equalizes,  cq  more  closely  approaches  unity  but  tQ  is 
increased  and,  for  the  assumption  (4.  2-10)  to  hold,  the  integration  time 
must  be  increased.  The  case  for  equalization,  then,  is  not  clear.  In  many 
practical  autocorrelation  systems  of  interest,  there  are,  then,  reasons  not 
to  equalize. 

Prior  to  the  formal  development  of  the  above  theory,  Mr.  Robert 
Bland  of  the  Department  of  Transportation's  Transportation  Systems  Center 
(TSC)  simulated  ranging  accuracy  as  a function  of  Aerosat  Filter  type  [only 
minimum  phase  filters  (Butterworth  and  Tchebyshev)  were  considered]. 
Although  the  effects  of  phase  and  amplitude  distortions  were  not  separated, 
Bland's  results  did  suggest  that,  in  accordance  with  the  above,  phase  dis- 
tortion, per  se,  had  little  effect  on  ranging  accuracy. 

(7) 

LeFande'  ' took  a frequency-domain  (as  opposed  to  the  present  time- 
domain)  approach  to  the  present  problem.  His  general  conclusions  are  the 
same  as  those  suggested  by  Bland's  work  and  some  preliminary  computations 
of  (4.  2-13)  (which  indicate  that  cq  < 1 for  a large  class  of  filters);  phase  dis- 
tortion is  probably  being  "tightly"  specified  to  an  extent  which  does  not  justify 
the  implied  expense  of  equalization. 

In  many  systems  of  practical  interest  (e.  g.  , HF /SSB  and  telephone 
modems),  the  luxury  of  integrating  over  many  signal  "chips"  is  not  afforded 
and  the  effects  of  differential  group  delay  are  more  complicated  and 
deleterious. 
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5.  0 DISCUSSION  AND  COMPARISON  WITH  PREVIOUS  WORK 

■ - .3, 

That  part  c_n  of  the  result  (2.  3-9)  is  the  same  as  that  of  Wheeler, 

but  that  author  considers  only  the  lowpass  case,  thus  obscuring  the  fact  that 

the  above  is  true  only  if  arithmetic  symmetry  prevails.  The  fact  that  arith- 

metric  asymmetry  leads  to  AM  to  PM  and  PM  to  AM  conversion  is  qualitatively 
(4). 

stated  by  Panter. 

The  work  in  Sections  3.  2 (in-band  amplitude  ripple),  and  that  in 

Sections  4.1  (differential  group  delay  of  the  Aerosat  wideband  channel)  and 

4.  2 (degradation  of  autocorrelation  signals)  are  believed  to  be  original.  The 

(6) 

latter  results  are  believed  to  be  more  general  than  those  stated  by  LeFande  , 

as  these  do  not  depend  on  a specific  signal  format,  but  only  on  the  method  of 

detection.  The  former  results  are  an  extension  of  the  work  of  Holmes  and 

(1) 

Keenan 

6. 0 CONCLUSIONS 

The  principal  conclusions  of  this  work,  stated  in  the  order  of  their 
perceived  importance,  are 

1)  For  any  linear  symmetric  bandpass  system  meeting  certain 
(practical)  constraints,  the  output  autocorrelation  signal  is 

related  to  the  input  autocorrelation  signal  by  a constant  which 
is  less  than  one  and  greater  than  zero.  This  constant  can  be 
computed  by  the  methods  given  in  this  work. 

2)  For  both  autocorrelation  and  "regular"  signals,  phase  non- 
linearity  is  best  specified  in  terms  of  t , as  opposed  to 

A t . 
gr 

3)  For  unequalized  system  transfer  functions,  the  effects  of  in- 
band  amplitude  distortion  on  phase  linearity  can  be  computed 
by  the  methods  of  Section  3.  2. 

This  (and  other)  work  suggests  that  "tight"  specification  of  phase  distortion  in 
systems  wherein  correlation  detection  is  employed  is  unnecessary. 
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Appendix  1 


Output  of  a Nonlinear  Phase  System 

The  purpose  here  is  to  derive  the  expressions  (2.  1-5)  and  (2.  1-7). 
Towards  this  end,  let 


-j  e(u;-wo) 


j2irn(u)  w3-cj2) 


tO->  < (O  < 10' 


where 


(A  1-  1) 


w3  " w2 


■jQ(w  _UJ0)  _ • 2r  n(ui  - to  -uj^) 

^ J 


(A  1-2) 


Then,  from  (2.  1-2),  that  part  of  the  frequency  response  due  to  system  phase 
can  be  written  as 


-j0(u>) 

e J ' = e 


■j[0o+to(w"wofl 


ao 


jnT(w-u>o) 


(A  1-3) 


~ ~ ft  u)  *C  to  oj  • 

w3  “ w2  r3  " f2  2 


The  impulse  response  of  the  phase  box  is  the  inverse  Fourier  transform  of 
the  frequency-domain  response  given  in  (Al-3): 


he(t*  ■ cn 


-j[e0-“0t0-nT"0]  ! f°  -jMt0+nr) 

2tF  / U'w>  e e dw* 


u>2  < to  < toj* 


(A  1-4) 


where  U(w)  is  as  given  in  Figure  2.  1-1  and  arises  due  to  the  finite  interval 
over  which  the  Fourier  expansion  of  ©(w  - u)q)  is  postulated  to  be 
applicable. 

Since  multiplication  in  the  frequency  domain  is  equivalent  to  convolu- 
tion in  the  time  domain,  (Al-4)  can  be  written  as 


-j  (O  ~u)  t -nTu  1 

h0(t)  = cn  e ° ° ° [u(t)  ® 6 (t  - tQ  - nr)],  (Al-5) 


where  u(t)  is  the  inverse  Fourier  transform  of  U(w)  and  use  has  been  made  of 


rcu 

+jw(t-to-nT) 


dw  = 6(t  - t - n^, 


(A  1-6) 


The  output  of  the  phase  box  in  response  to  an  input  a(t)  e^w*  is  just  the  con- 
volution of  the  input  with  the  impulse  response  of  the  phase  box; 


yQ(t)  = x(t)  ® hQ(t) 


= a(t)  ejwt  ® ^ cn 


-if©  -w  t -nTw  "1 
JLuo  o o o-l 


u(t)  ® 6(t  - tQ  - nt) 


■IN 


j[wt-eo+(Wo-w)(to+nT)] 


a(t  - tQ  - nT)  ® u(t) 


(A  1-6) 


where  the  order  of  convolutions  has  been  changed  and  use  has  been  made  of 
the  property  f(t)  ® 6(a)  = f(a). 
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The  output,  y(t),  of  the  entire  system  (including  the  magnitude  box) 
is  just  (Al-6)  convolved  with  the  impulse  response  of  the  magnitude  box, 


h (t): 
m 


y(t) 


CO 

E 


Ln=  -co 


c e 
n 


j[wt_eo+(wo~w^to+nT^l 


a(t  - t - nr) 


® u(t) 


= ye(t)  ® u(t)  ® hm(t) 

= ye(t)  ® |H(w)|u(w) 


(A  1-7) 


where 


ye(t) 


CO 

-£ 


c e 
n 


j[wt-eo+(Wo-w)(to+nT)] 


a(t  - tQ  - nT) 


(A  1-8) 


n=  -a) 


The  results  (Al-6)  and  (Al-7)  are  those  stated  in  (2.  1-5)  and  (2.  1-7) 
respectively. 
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Appendix  2 


Group  Delay  Due  To  In-Band  Amplitude  Distortion 

The  initial  purpose  here  is  to  derive  the  results  (3.  2-2)  and  (3.  2-5). 
Towards  this  end,  it  is  noted  that,  for  a minimum-phase  system,  the  phase 
is  related  to  the  amplitude  characteristics  by  the  Hilbert  transform:^ 


6(u>)  = ^ I In 


(5) 


w - x d L I H(x)l 


+ x dx 


, I H(x)l  , 

ln  I ■ “ 


2 < W < W^»  (A2-  1) 


Using  Leibnitz'  rule'  , the  group  delay  is 


t (w) 

gr'  ' 


de(w) 

du> 


/ ( 

l I 7Z  !ln 


Ul  - X 


W + X 


dx 


ln 


Ih^1 


|H(u>2) 


'dx, 


(A2-2) 


u>2  < w ui^. 


For  the  arbitrary  amplitude  distortion  cited  in  (3.  2-1), 


as 


[ln  F^i] = £ 

L J n=l 


7rnb 


cos  nn 


Wo  " W- 


|~x  ~ w2l 

[W3  - w2J’ 


(A2-3) 


and  (A2-2)  becomes* 


nb 




w,  - , 


n=l 


w - x 


cos  nir 


fX  " “2! 

l“3  ‘ w2j 


dx 


(A2-4) 


w + x 


cos  nn 


I x ~ w2  1 

Lw3  ‘ w2j 


dx 


Wo  < W < W- 


* Removal  of  the  magnitude  sings  in  lu>-xl  in  (A2-2)  has  been  effected  by 

I U>3)  S° 
ting  integral 


breaking  the  interval  of  integration,  ( ou  , w ),  into  (w  , u>)  and  (u>  a)  ) so  as 

• j ! ^ _ * a.1  a * ~ ml  1 A*  J • 

resultii 


to  maintain  111 -x  or  x - id  orx  • u)  positive  quantities.  The 
is  given  in  A2-4). 

31 

I --  — I ■ 


Now,  both  the  variable  x and  the  band  of  interest  (in  w)  is  confined  to 
(u)2»  w-j).  If  this  interval  is  small  in  comparison  to  u>q  (=  (w^  + the 

magnitude  of  the  second  integral  in  (A2.  4)  is  on  the  order  of  - w^)  /2ojq 
times  the  magnitude  of  the  first  integral.  Hereafter,  it  is  assumed  that 

2u>o  » (w^  - (A 2 


and  the  second  integral  in  (A2-4)  is  neglected. 

With  the  change  of  variable  £ = oj  - x in  the  first  integral  in  (A2-4), 


00 

E 


nb 


n 


t (oj)  a - 

gr'  u>3  - co2 

n=l 


/•OJ-OO-3 

/ 3 1 

_e  u “ w2 

1 —t—  cos  ntT 

2 + 

! r 

w3  ‘ w2  w3  w2 

0J«>  < OJ  < OJ^, 


2u>^  » (oj^  - ^2)' 


(A2- 


Using  cos  (x  + y)  - cos  x cos  y - sin  x sin  y,  (A2-6)  can  be  written  as 


With  the  further  change  of  variable 


(A2-7)  can  be  written  as 


t (w)  = 
gr'  ' 


^ nb 

Y — 

/ i Wo  - <J 


n=l 


3 " w2 


cos  X. 

^ dX 


sinX 


dX 


(A2-9) 

By  letting  w = wQ  + Aw,  expanding  the  cos  and  sin  functions  (using  well-known 
expressions  for  cos(x  + y),  sin(x  + y))  preceding  the  first  and  second  integrals 
in  (A2-9),  respectively,  and  making  use  of,  e.  g.  , 


- cos  mr/2  = 


.1 


, n odd 


|(-d 


!♦* 


(A2-10) 


n even 


it  can  be  shown  from  (A2-9)  that 


n , 

2+1k 
co  ( - 1 ) nb 


’ -a+bAw 


tgrK+Aw)  - 
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nirAw 


/ , \ \ sin 

( 3 - W2>  ( w3  “ w2  , 
n-  1,  3*  • • ' •'a-bAui 


cos  X 
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•'a+bAw 


sinX 
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(A2-11) 

(cont’d) 
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§+l 

» (-l)2nb  ( r-a+ba“  , 

+ > 'cos I — r — dX 

4-i  w3  - w2  w3  - w2  I , . X 

n=2,  4*  • * V •'a+bAu; 


where 


nir  , 
a = -x— , and 


(A2-12) 


w3  " w2 


No  approximation  other  them  (A2-5)  has  been  made  in  deriving  the  ~bove 


result. 


Reduction  of  (A2-11)  to  (3.  2 - 2)  is  straight  forward. 


. nirAw 
+ sin 

w3  " w2 


- a+bAw 


a+bAu) 


sin  X , . I 

dxj' 


(A2-11  cont'd) 
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In  the  arithmetically  symmetric  case, 

b = o n even  and  the  terms  following  the  second  £ are  zero.  Then,  it 
n ’ 

follows  from  the  above  that 


t (w  - Aw)  = t (w  + Aw) 
gr  o gr  o • 


(A2-25) 


This  is  the  result  stated  in  (2.3-5).  Finally,  it  should  be  noted  that  in  pres- 
ent development,  the  interval  (w2,  w3)  has  not  been  restricted  to  the  "bandpass" 
or  any  other  segment  of  the  response  of  a system.  Therefore,  the  afore- 
mentioned results  are  generally  applicable. 


Appendix  3 


Fourier  Expansion  of  In-Band  Amplitude  Distortion 


If  the  in-band  amplitude  distortion 


F(w  - w_)  = In 


[JM1 
[ih(W2)|  J 


(A3-  1) 


has  a period  twice  the  interval  and  is  odd  about  (he.,  T(uj-uj^)  = 

- F(oj2  - oj)), 


F(  OJ  - o,2) 


CO 

= > b sin  nu 

Z_J  n 


W3  - w2 


(A3-  2) 


where 


! I 


b = -57 — i r 

n 2(w3  -w2) 


F( oj  - oj  2)  sin  nTr 


2^2"w3 


(=Pt!  -• 


°3  r 2 r 3 

[•••]  = I [•••]  + / [*  • *] 


2 u)  ^ -co  3 ^ 2oj^-cj3 


where  • • • indicates  the  kernel  of  (A 3- 3).  Since  F(oj  - u>2)  is  odd* 


[•••]=  / [•••], 


2u)2~<jJ3 


(A3-3) 


(A3-4) 


(A3-5) 
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so  that 


n 


U)_  _ GJ  « 

3 2 


F(  w - u>0)  sin  Tin  — 

C u> 


3 "w2 


(A3-6 


This  is  given  in  (3.  2-la). 

Since  (A3-2)  is  an  adequate  representation  of  F(w  - w^)  in  the  inter- 
val (2w-,  - it  is  an  adequate  representation  of  F(w-u^)  in  the  lesser 

interval  (u^.u^).  Thus,  (A3-2)  is  used  to  represent  F(u-w^)  and  is  applied 
in  the  bandpass  interval  (w2,co^)  only. 


Appendix  4 


Noise  Bandwidth  of  an  Idealized  Filter 


< 


I * 

Ri 

1 ! 

I 4 


The  purpose  here  is  to  prove  the  result  (4.  1-1).  This  is  applicable 
to  a system  with  the  characteristics  shown  in  Figure  1.0-1  with  a flat  in-band 
response  and 


n 


(A4-  1) 


In  the  region  w ^ < u>  < w 


In  |H(w)| 


(A4-2) 


or 


a[w-w  . /H] 

|H(«)|  = |H(w1)J  e 


(A4-3) 


where 


(A4-4) 


Similarly,  in  the  region  w^<  w 


I H(w)| 


- a [u>  - 

e 


/Hi 


(A4-  5) 


1 

i ' 


vi 


Because  a is  chosen  to  be  large  in  the  numerical  example,  the  contribution 
to  the  noise  bandwidth  in  the  regions  w > and  u>  <■  is  neglected;  this 
amounts  to  neglecting  contributions  to  the  noise  bandwidth  which  lie  outside  of  the 
60dB  points  of  a practical  apprximation  of  the  idealized  magnitude  charac- 
teristic considered. 
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With  the  above  assumption,  the  noise  bandwidth  is 


^ry  |H(u>)|2  du>  = Uy  / 

2>|ZX  |h<“2>|2X1 


H(wi)|2  fw 2 2a[w-w1  /n] 


/*w3  fw4  -Zatw-uij/n] 

+ ( dw  +1  e dw 


By  carrying  out  the  integrations  above,  this  is  readily  shown  to  be 


N - (w.j  w2)  +f, 


|H(Wl) 

|H(«2) 


This  is  the  result  given  in  (4.  1-1). 


40 


